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Abstract 

We  study  the  heavy  traffic  regime  of  a  discrete-time  queue  driven  by  cor¬ 
related  inputs,  namely  the  M\GI\oo  input  processes  of  Cox.  We  distinguish 
between  M|GI|oo  processes  with  short-  and  long-range  dependence,  identi¬ 
fying  for  each  case  the  appropriate  heavy  traffic  scaling  that  results  in  non¬ 
degenerate  limits.  As  expected,  the  limits  we  obtain  for  short-range  dependent 
inputs  involve  the  standard  Brownian  motion.  Of  particular  interest  are  the 
conclusions  for  the  long-range  dependent  case:  The  normalized  queue  length 
can  be  expressed  as  a  function  not  of  a  fractional  Brownian  motion,  but  of  an 
a-stable,  1/a  self-similar  independent  increments  Levy  process.  The  resulting 
buffer  asymptotics  in  heavy  traffic  display  a  hyperbolic  decay,  of  power  1  -  a. 
Thus,  M\GI\oo  processes  already  demonstrate  that,  within  long-range  depen¬ 
dence,  fractional  Brownian  motion  does  not  necessarily  assume  the  ubiquitous 
role  that  standard  Brownian  motion  plays  in  the  short-range  dependence  setup. 


1  Introduction 

The  apparent  presence  of  long-range  dependence  and  self-similarity  in  network  traf¬ 
fic  has  been  suggested  by  several  traffic  measurement  studies  (e.g.,  WAN  [24],  Ether- 
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net  [17]  and  VBR  video  [3]),  and  points  to  the  need  of  revisiting  various  performance 
analysis  and  design  issues.  Indeed,  as  some  recent  experimental  work  [10]  already 
suggests,  long-range  dependence  has  a  tangible  and  adverse  effect  on  queueing  mea¬ 
sures  such  as  buffer  overflow  probabilities,  and  its  presence  thus  cannot  be  overlooked 
or  underestimated. 

Roughly  speaking,  long-range  dependence  amounts  to  correlations  in  the  packet 
stream  spanning  multiple  time  scales,  which  are  individually  rather  small  but  de¬ 
cay  so  slowly  as  to  be  non-summable.  This  is  expected  to  affect  performance  in  a 
manner  drastically  different  from  that  predicted  by  (traditional)  summable  corre¬ 
lation  structures  which  typically  arise  in  Markovian  models.  In  order  to  gam  some 
understanding  into  long-range  dependence,  and  its  impact  on  queueing,  various 
traffic  models  with  non-summable  correlation  patterns  have  been  proposed  m  the 
literature.  These  models  include,  among  others,  fractional  Gaussian  noise  inputs 
[1],  fractional  Brownian  motion  [21]  and  on/off  sources  with  subexponential  inter¬ 
renewal  periods  [7,  14]  .  In  all  these  cases,  the  corresponding  buffer  asymptotics  do 
not  display  the  exponential  tails  typically  associated  with  short-range  Markovian 

models. 

In  this  paper,  we  consider  the  class  of  discrete-time  M\GI\oo  input  processes.  An 
M\GI\oo  input  process  is  understood  as  the  busy  server  process  of  a  discrete  time  in¬ 
finite  server  system  fed  by  a  discrete-time  Poisson  process  of  rate  A  (customers/slot) 
and  with  generic  service  time  a.  Such  M|GJ|oo  processes  can  account  both  for 
short  and  long-range  dependent  behaviors,  with  the  correlation  patterns  controlled 
through  cr  [Proposition  2.1],  Furthermore,  asymptotic  self-similarity  arises  when  cr 
is  Pareto-like,  i.e.,  has  a  regularly  varying  tail  of  the  form  (4.1).  M\GI\oo  processes 
have  already  been  used  by  Paxson  and  Floyd  to  successfully  model  WAN  traffic  [24]. 
However,  the  relevance  of  M\GI\oo  input  processes  to  network  traffic  modelling  is 
perhaps  best  explained  through  the  fact  they  naturally  arise  as  the  aggregate  limit 

of  a  large  number  of  on-off  sources  [18]. 

As  discussed  in  [9,  19,  22,  23],  M\GI\oo  processes  induce  a  wide  variety  of 
asymptotic  behaviors  for  the  buffer  probabilities  at  a  multiplexer  with  constant 
release  rate.  In  particular,  when  a  has  a  regularly  varying  tail  -  the  M\GI\oo 
process  is  now  asymptotically  self-similar,  the  buffer  asymptotics  are  hyperbolic  in 
nature,  in  stark  contrast  with  the  Weibullian  tails  induced  by  fractional  Gaussian 
noise  (or  fractional  Brownian  motion)  [21].  Here,  we  further  explore  this  discrepancy 
in  the  heavy  traffic  regime  [12],  In  fact,  the  motivation  for  going  to  heavy  traffic  is 
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manifold: 

First,  under  short-range  dependence  heavy  traffic  analysis  has  offered  useful 
characterization  of  queueing  networks  in  terms  of  functionals  on  Brownian  motion 
[12],  and  it  is  therefore  of  theoretical  interest  to  extend  the  analysis  to  the  long-range 
dependence  setup.  Such  an  extension  would  possibly  help  answer  the  question  as  to 
whether  fractional  Brownian  motion,  the  long-range  dependent  analog  of  standard 
Brownian  motion,  does  play  a  similar  central  role  in  the  modeling  of  long-range 

dependent  traffic. 

Next,  heavy  traffic  information  constitutes  a  key  component  of  the  light  traffic 
interpolation  technique  originally  developed  by  Simon  and  Reiman  in  a  series  of 
papers  [25,  26,  27]  in  order  to  estimate  performance  metrics  such  as  response  times 
in  queueing  systems  driven  by  Poisson-like  inputs.  As  we  contemplate  the  possible 
use  of  this  technique  with  M|G/|oo  input  streams,  we  are  naturally  led  to  seek  a 
complete  classification  of  the  heavy  traffic  limits  which  arise  when  considering  such 

inputs. 

Finally,  a  heavy  traffic  analysis  of  Af  |G/|oo  processes  might  help  elucidate  the 
noted  difference  in  buffer  asymptotics  between  M\GI\oo  and  fractional  Gaussian 
noise  inputs.  Unless  this  is  due  to  some  fundamental  structural  property,  both 
models  are  expected  to  have  a  heavy  traffic  characterization  in  terms  of  fractional 
Brownian  motion,  in  very  much  the  same  manner  that  different  short-range  depen¬ 
dent  models  eventually  collapse  to  a  single  description  involving  Brownian  motion 
[12].  Otherwise,  despite  asymptotically  identical  correlation  patterns,  the  differ¬ 
ences  would  carry  over  and  manifest  themselves  even  more  clearly  m  the  heavy 

traffic  regime. 

The  results  presented  here  confirm  the  latter  possibility:  Under  short-range 
dependence,  the  class  of  Af|G/|oo  inputs  is  found  to  belong  to  the  domain  of  at¬ 
traction  of  the  standard  Brownian  motion,  as  expected.  More  significantly,  we  show 
that  under  long-range  dependence,  with  a  belonging  to  the  domain  of  attraction 
of  a  non-normal  stable  law,  the  M\GI\oo  process  is  not  attracted  to  a  fractional 
Brownian  motion,  but  instead  to  a  non- Gaussian,  a-stable  Levy  motion  which  is 
1/a  self-similar.  As  a  consequence,  the  distribution  of  the  heavy  traffic  queue  length 
does  not  display  a  Weibullian,  but  a  Pareto  tail,  with  power  1  -  a  [Theorem  5.3]. 
These  results  underscore  the  fundamentally  different  nature  of  the  long-range  de¬ 
pendent  M\GI\oo  process  (when  compared  to  fractional  Gaussian  noise),  and  also 
point  to  the  fact  that  fractional  Brownian  motion  does  not  necessarily  play  for  long- 
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range  dependence  the  same  key  role  that  standard  Brownian  motion  assumes  under 
short-range  dependence.  Within  long-range  dependence,  there  seems  to  be  a  choice 
for  distinct  modeling  possibilities,  and  it  is  not  at  all  difficult  to  find  rather  simple, 
potentially  useful  traffic  models  that  are  attracted  to  non-Gaussian  limits. 

The  basic  idea  behind  the  proof  of  these  results  is  a  “convergence  together”  argu¬ 
ment  which  allows  us  to  identify  processes  with  well-known  heavy  traffic  behavior, 
under  both  short-  and  long-range  dependence.  This  is  accomplished  chiefly  by 
combining  standard  results  on  stable  random  variables  and  their  domain  of  attrac¬ 
tion  [11],  with  a  general  functional  convergence  result  for  processes  with  stationary 
independent  increments  due  to  Skorokhod  [29].  We  point  out  that,  even  in  the 
short-range  dependent  case,  convergence  to  Brownian  motion  does  not  appear  to 
follow  from  standard  results  for  stationary  processes  [4,  Thm.  20.1,  p.  174],  as  it 
is  not  obvious  that  the  M\GI\oo  busy  server  process  satisfies  the  required  mixing 
property.  Note  however  that,  as  shown  in  [22],  the  M\GI\oo  busy  server  process 
is  strongly  positively  correlated  -  it  is  an  associated  process.  Because  of  this  prop¬ 
erty,  it  is  then  possible  under  short-range  dependence  to  develop  an  alternative 
approach  [32],  similar  to  that  used  by  Newman  and  Wright  in  [20]  in  establishing 
the  Invariance  Principle  for  associated  random  variables. 

Related  work  on  heavy  traffic  queueing  analysis  under  long-range  dependence 
appears  to  have  been  initiated  by  Norros  [21],  where  the  presence  of  fractional  Brow¬ 
nian  motion  is  postulated.  This  line  of  inquiry  is  further  pursued  in  [30],  while  in 
[7]  Brichet  et  al.  show  how  fractional  Brownian  motion  can  arise  from  a  Gaus¬ 
sian  superposition  scheme  of  infinitely  many  on/ off  sources  with  heavy  tailed  on/off 
periods.  In  view  of  the  fact  that  M\GI\oo  processes  arise  from  a  different  super¬ 
position  scheme  of  infinitely  many  on/off  sources  [18],  it  is  not  too  surprising  that 
these  lead  to  a  different  heavy  traffic  limit,  involving  Levy  motions.  Heavy  traffic 
results  similar  and  related  to  the  ones  given  here  have  been  reported  in  [15],  where 
only  convergence  of  finite  dimensional  distributions  is  announced.  The  conclusions 
discussed  here  were  obtained  independently,  and  were  summarized  in  the  conference 
paper  [31]. 

The  remainder  of  the  paper  is  organized  as  follows:  The  class  of  M\GI\oo  input 
processes,  along  with  the  discrete-time  queueing  setup,  are  introduced  in  Section 
2.  We  explain  how  the  queue  is  driven  to  heavy  traffic  in  Section  3.  The  main 
heavy  traffic  results  are  then  stated  in  Section  4,  while  Section  5  discusses  their 
consequences  on  the  queue  length  asymptotics.  An  outline  of  the  proofs  is  presented 
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in  Section  6.  The  proofs  of  the  main  results  follow  in  Section  7,  while  the  arguments 
behind  the  “convergence  together”  are  discussed  in  Section  8.  In  the  appen  ix 
Sections  9  and  10,  we  have  summarised  several  technical  facts  concerning  unc  ions 

A  few  words  about  the  notation  used  here.  All  rvs  are  defined  on  some  probabil¬ 
ity  triple  (fi,  T,  P),  with  E  denoting  the  corresponding  expectation  operator.  We  use 
=>r  to  denote  weak  convergence  [4],  and  Ar  to  denote  convergence  m  probability 
(with  r  going  to  infinity).  We  write  f(x)  ~  g{x)  (x  ->•  oo)  when  Jim  f(x)/g(x)  . 

Equality  in  distribution  is  denoted  by  =at,  and  inequality  in  the  strong  stochas  ic 

order  is  denoted  by 

2  A  multiplexer  driven  by  M\GI\oo  input  processes 

We  begin  by  collecting  some  pertinent  facts  concerning  the  busy  server  process  of 
a  discrete-time  M|G/|oo  system;  the  reader  is  referred  to  [8,  22]  for  proo  s  an 

additional  material  on  this  class  of  processes.  ,  , 

Consider  a  system  with  infinitely  many  servers:  During  time  slot  [n, n  +  ),  n 

/3n+1  new  customers  arrive  into  the  system.  Customer  j,  j  -  1 . Pn+u 

is  presented  to  its  own  server  and  begins  service  by  the  start  of  slot  [n  +  l,n  +  2); 
its  service  tune  has  duration  a„+lj  (in  number  of  slots).  Let  6,  denote  the  number 
of  busy  servers,  or  equivalently,  the  number  of  customers  still  present  m  the  system 
at  the  beginning  of  time  slot  [n,  n  +  D-  If  initially  (he.,  at  time  n  =  0)  there  were 
already  b  busy  servers,  we  denote  by  aj  the  residual  service  duration  (m  time  slots) 

for  the  jth  busy  server,  j  =  l . 6-  The  busy  server  process  {bn,  n-  0, 1, .  •  is 

what  we  refer  to  as  the  M\GI\oo  input  process.  .  _ 

Throughout,  the  W-valued  rvs  b,  {/?„+!,  n  =  0, 1, . .  n=  1,  ,  •  •  •  5  3  ~ 

12  1  and  {a  j  =  1,2,...}  satisfy  the  following  assumptions:  (i)  These  rvs 

are' mutually  independent;  (11)  The  rvs  "  =  0,1,...}  are  i.U.  Poisson  rvs 

with  parameter  A  >  0;  (ill)  The  rvs  {<r„,„  n  =  1, . . .  1  j  =  U,  ■  ■  •>  «  «•••*  "th 
common  pmf  G  on  {1,2, . .  .}■  Let  cr  be  a  generic  K^valued  rv  distributedaccordmg 

to  the  pmf  G,  assume  throughout  that  E  [a]  <  oo;  (iv)  The  rvs  {ctj,  ^  " 

are  *  *  d.  W-valued  rvs  distributed  according  to  the  equilibrium  pmf  G  associated 
with  G,  i.e.,  if  a  denotes  a  generic  W-valued  rv  distributed  according  to  the  pmf 

then  P  [o>n\ 

E  [cr] 


p  [a  =  n]  = 


n  =  1,2, . . . 


(2.1) 
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The  tail  of  the  equilibrium  pmf  G  directly  controls  the  correlations  of  the  se¬ 
quence  { bn ,  n  =  0, 1, .  ■ 

Proposition  2.1  If  b  is  taken  to  be  a  Poisson  rv  with  parameter  AE[o),  then  the 
busy  server  process  {&„,  n  =  0, 1, . . .}  is  a  (strictly)  stationary  ergodic  process  mth 

the  following  properties: 

1.  For  each  n  =  0, 1, . . the  rv  bn  is  a  Poisson  rv  with  parameter  AE  [o J; 

2.  Its  covariance  function  is  given  by 

cov(f>„+j,M  =  AE  [[<7  -  j]+]  =  AE  M  P  [if  >  ]} ,  ".j  =0,1,... 

3.  Its  index  of  dispersion  of  counts  (IDC)  is  given  by 

oo  00  A 

IDC  =  y;  cov(6„+j,  bn)  =  AE  [a]  £  p  [5  >  j]  =  ^  [°{p  +  ^  > 
j=o  i=° 

and  the  process  is  short-range  dependent  (i.e.,  IDC  finite)  if  and  only  if  E  [a2]  is 
finite. 

In  short,  a  stationary  M|G/|oo  input  process  is  fully  characterized  by  the  pair 
(A,  a),  and  displays  time  dependencies  which  are  determined  by  the  tail  of  a.  We  now 
offer  such  a  stationary  M|GI|oo  input  process  {bn,  n  =  1,2,... , }  to  a  multiplexer 
which  we  model  as  a  discrete-time  single  server  queue  with  infinite  buffer  capacity, 
operating  at  a  constant  rate  and  in  a  first-come  first-served  manner.  Let  qn  denote 
the  number  of  cells  remaining  in  the  buffer  by  the  end  of  slot  [n  -  1,  n),  and  let  bn+ 1 
denote  the  number  of  new  cells  which  arrive  at  the  start  of  time  slot  K  n  + 1) .  If  the 
multiplexer  output  link  can  transmit  c  cells/slot,  then  the  buffer  content  sequence 
n  =  0, 1,  ■  •  •}  evolves  according  to  the  Lindley  recursion 

q0  =  0;  qn+i  =  [<ln  +  bn+i  -  c}+ ,  n  =  0,1,...  (2-2) 

By  Part  1  of  Proposition  2.1,  the  average  input  rate  to  the  multiplexer  is  simply 
E  [bn\  =  AE  [ct],  and  it  can  be  shown  that  the  system  is  stable  if  AE  [a]  <  c,  in  which 

case  qn  =>n  <7oo  for  some  IN  -valued  rv  r/oc- 

The  output  to  the  Lindley  recursion  (2.2)  admits  an  equivalent  representation, 

which  is  useful  for  establishing  heavy  traffic  limit  theorems:  We  can  write 
qo  =  0;  qn  —  sn  ~  nc  —  inf  (sj  -  jc,  j  =  0, 1,  •  ■  •  ,  n) ,  n  -  1, 2, . . . 

upon  defining  the  partial  sums  {sn,  n  =  0, 1,  •  •  ■}  by 

s0  =  0;  n  =  l,2,...  (2-3) 

i- 1 
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3  The  heavy  traffic  regime 

We  seek  to  understand  the  behavior  of  the  (stable)  queue  under  the  assumption  that 
it  is  almost  fully  utilized,  i.e.,  AE  [cr] ,  though  less  than  the  release  rate  c,  is  very 
close  to  c.  This  typically  involves  obtaining  limiting  expressions  of  properly  rescaled 
quantities  of  interest,  as  the  traffic  intensity  p  =  AE  [er]  /c  tends  towards  its  critical 
value  1.  Here,  the  quantity  of  interest  is  the  steady-state  queue  size  (foe-  A  natural 
setup  to  investigate  this  problem  consists  of  embedding  the  discrete-time  queue  with 
release  rate  c  driven  by  an  M\GI\oo  input  process  (A,er)  into  a  parametric  family 
of  like  queueing  systems,  indexed  by  an  integer  parameter,  say  r.  More  precisely, 
for  each  r  =  1,  2, ...  we  take  the  rth  system  to  be  a  discrete-time  queue  with  release 
rate  c  driven  by  an  M\GI\oo  input  process  {&(),  n  =  0, 1, . . .}  characterized  by  the 
pair  (Ar,  cr).  The  corresponding  queue  size  sequence  { qrn ,  n  =  0, 1, . . .}  also  obeys 
the  Lindley  recursion  (2.2),  and  can  be  represented  as 

qrn  =  srn  -  nc-  inf  (srj  -  jc,  j  =  0, 1,. ..  ,n)  ,  n  =  l,2,...  (3.1) 

where  {s(),  n  —  1,2, . . .}  is  the  sequence  of  partial  sums  (2.3)  associated  with 
{brn,  n  =  1,2,...}.  We  take  ArE  [<r]  <  c  for  all  r  =  1, 2, . . with 

lim  Ar  =  c/E  [cr] .  (3.2) 

r— >oo 

Thus,  each  one  of  these  systems  is  stable  with  lim  pr  =  1,  thereby  capturing  the 
notion  that  “the  system  is  driven  to  heavy  traffic.”  We  seek  a  scaling  sequence 
{(r,  r  =  1,2,...}  such  that  the  convergence  in  distribution 

f  =>rQ  (3.3) 

Cr 

takes  place  to  some  IR-valued  rv  Q. 

Unfortunately,  this  heavy  traffic  program  cannot  be  carried  out  in  this  form  as 
exact  expressions  are  unavailable  for  the  distribution  of  q ^  owing  to  the  correlations 
present  in  the  M\GI\oo  input  process,  and  we  need  to  resort  to  the  following  indirect 
approach  where  the  buffer  content  is  rescaled  in  both  the  time  and  state  space 
variables.  Thus,  for  each  r  =  1,2,...,  we  define  the  IR-valued  continuous-time 
processes  {Sr(t),  t  >  0}  and  {Qr{t),  t  >  0}  by 

T 

S’W  =  1-  (4„  -  E  [4,])  and  =  t>  0, 
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and  we  introduce  the  function  'f  :  IR.-^  — >  IR  given  by 

7r«)  =  Y  (Me  -  E  fSflrJ)  =  M  (c  -  ArE  [tr]) ,  t  >  0. 

Sr  >r 

We  note  that  (3.3)  can  informally  be  stated  as 

lim  limQr(f)  =  Q  (3-4) 

1 - >OOt— >00 

with  limits  understood  in  the  sense  of  weak  convergence.  The  approach  to  heavy 
traffic  followed  here  is  to  interchange  the  order  of  these  limits,  i.e.,  to  evaluate 

lim  lim  Qr(t)  (3-5) 

«->oo  r— HX> 

which  corresponds  to  first  taking  r  — >  oo,  and  then  letting  t  — >  oo.  Assuming  that 
the  limits  can  be  taken  in  that  order,  we  are  then  left  with  the  task  of  showing  that 


lim  lim  Qr(t)  =  Q  =  lim  lim  Qr{t). 

T— >00  t— ¥00  t— >oor— HX> 


(3.6) 


In  this  paper  we  concentrate  only  on  establishing  the  first  step  (3.5),  and  it  is 
well  known  [12,  33]  that  the  theory  weak  convergence  on  function  spaces  provides 
a  natural  framework  for  doing  so.  To  that  end,  we  pause  briefly  to  introduce  the 
needed  notation,  as  well  as  to  highlight  several  points  from  the  theory  of  weak 
convergence  of  processes;  this  material  is  drawn  mostly  from  [4,  pp.  150  153]  to 
which  the  reader  is  referred  for  additional  information: 

For  each  T  >  0,  let  £>[0,T]  denote  the  space  of  mappings  [0,T]  -4  IR  which 
are  right-continuous  with  left  limits;  the  space  D[0,T]  can  be  equipped  with  either 
the  uniform  topology  or  the  standard  Skorokhod  topology  [4,  p.  111].  As  in  [4,  p. 
150],  a  concept  prefixed  with  U  (resp.  S )  refers  to  the  uniform  (resp.  Skorokhod) 
topology.  For  probability  measures  defined  on  the  collection  of  [7-Borel  (resp.  S- 
Borel)  sets  on  D[0,T],  we  refer  to  weak  convergence  in  the  sense  of  the  uniform 
(resp.  Skorokhod)  topology  by  C/-weak  (resp.  S-weak)  convergence,  and  we  write 
=>r  (resp.  ===>r)  (with  the  understanding  that  r  goes  to  infinity).  For  probability 
measures  defined  on  the  collection  of  (7-Borel  sets,  [/-weak  convergence  implies  S- 
weak  convergence  but  the  converse  is  false.  This  implication  will  be  used  repeatedly 
in  various  technical  arguments  [Sections  6  and  8]. 

Finally,  let  D[0,  oo)  denote  the  space  of  mappings  R+  ->  IR  which  are  right- 
continuous  with  left  limits.  In  this  paper,  we  present  results  on  the  5-weak  con¬ 
vergence  of  the  restrictions  to  finite  intervals  of  sequences  of  IR-valued  processes 
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with  sample  paths  in  Z?[0,oo).  More  precisely,  consider  the  sequence  of  It-  valued 
processes  {XT(t),  t  >  0},  r  =  1, 2, . . with  sample  paths  in  D[ 0,  oo).  Whenever,  for 
each  T  >  0,  we  have  the  S-weak  convergence 

{Xr{t),  0  <  t  <  T}  =Ur  {X(t),  0  <  t  <  T}  mD[0,T]  (3.7) 

for  some  IR-valued  process  {A(t),  t>  0}  with  sample  paths  in  I>[0,  oo),  we  simplify 
the  notation  by  writing 

{Xr(t),  t  >  0}=>r{X(t),  t  >  0}.  (3.8) 

Now,  noting  that  (3.1)  can  be  rewritten  as 

Qr(t)  =  Sr(t)-Y(t)~  inf  (Sr(x)-Y(x)),  i>0,  (3.9) 

0<x<t 

and  recalling  the  continuous  mapping  theorem  [4,  Thm.  5.1,  p.  30],  we  conclude 
that  the  first  limit  in  (3.5)  requires  at  the  very  least  identifying  a  scaling  sequence 
{(r,  r  —  1,  2, . . .}  that  ensures  the  convergence 

{Sr{t),  t>0}=>r  { S(t ),  t  >  0}  (3.10) 

for  some  non-trivial  limiting  process  {S(t),  t  >  0}. 


4  The  main  heavy  traffic  results 

As  will  become  apparent  shortly,  the  choice  of  the  scaling  sequence  {Cr,  r  =  1,  2, . . .} 
and  the  characterization  of  the  limiting  process  {S(t),  t  >  0}  entering  (3.10)  both 
depend  on  the  distribution  of  the  rv  a  which  controls  the  correlations  in  the  input 
cell  stream.  It  is  nevertheless  easy  to  see  that  in  order  to  avoid  collecting  only  a  law 
of  large  numbers,  any  candidate  scaling  sequence  {(r,  r  =  1,2,...}  should  obey  the 
following  necessary  condition: 


Condition  (A)  The  scaling  sequence  {Cr,  r  =  1,2, . . .}  satisfies 

Cr 

lim  Cr  =  +°o  with  lim  —  =  0. 

1 — too  r— Kxj  f 


The  heavy  traffic  assumption  below  refines  (3.2),  and  guarantees  that,  as  r  — >  oo, 
the  family  of  queueing  systems  described  by  (3.9)  gradually  approaches  instability 
at  the  appropriate  speed: 
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Assumption  (A)  The  scaling  sequence  {Cr,  r  =  1, 2, . . .}  satisfies 

lim  (ArE  [crl  -  c)  - 7  or  equivalenty ,  ArE  [c a]  =  c  -  ~(nf  +  o(l)) 

r-y  00  C,r  T 

for  some  7  >  0. 


Condition  (A)  and  Assumption  (A)  are  enforced  throughout.  It  is  worth  pointing 
out  that  the  scaling  sequence  {Cr,  r  =  1,2, . . .}  is  essentially  unique,  i.e.,  any  other 
scaling  sequence  {£,  r  =  1,2,...}  yielding  a  non-degenerate  limit  in  (3.10)  must 
satisfy  lim  C/Cr  =  C  for  some  finite  constant  C  >  0. 

We  begin  with  the  case  where  the  M\GI\oo  process  is  short-range  dependent 
and  let  {B(t),  t>  0}  denote  a  standard  Brownian  motion. 

Theorem  4.1  (Short-range  dependence)  If  E  [cr2]  <  00,  then  with  Cr  =  sfr,  r  = 
1,2,...,  it  holds  that 

(S’-(t),  t>  0}  (>0(. 

The  remaining  results  are  obtained  under  the  additional  assumption  that  the 
tail  of  cr  is  regularly  varying  of  order  a  (1  <  a  <  2),  i.e.,  of  the  form 


P  [cr  >  n]  =  n  ah{n),  n  =  0,1,... 


(4.1) 


for  some  slowly  varying  function  h  :  1R_)_  — >  1R+  [11],  in  which  case  the  mean  E  [a] 
is  finite.  Of  particular  interest  for  the  forthcoming  discussion  is  the  realization  that 
the  truncated  second  moment  of  cr  is  (2  -  a)-regularly  varying.  Writing 

a 


^a(^)  —  N 


h(x )  if  1  <  a  <  2 

2  —  a 

h{r) 

^2-j  r 


if  a  =  2 


(4.2) 


r= 1 


for  all  x  >  0,  we  can  show  via  Proposition  9.2  that  the  function  la  :  IR+  — >  1R+  is 
slowly  varying  and,  whenever  E  [cr2]  =  00, 


E  [l  [cr  <  n]a2j  ~  n2  Q/a(n),  (n  ->  00).  (4.3) 

The  details  of  the  proof  of  this  asymptotic  equivalence  are  identical  to  those  of  (7.10) 
and  (7.11). 

The  next  proposition  handles  the  boundary  value  a  —  2,  which  represents  a 
hybrid  case  between  short  and  long-range  dependence. 


10 


Theorem  4.2  Assume  a  =  2  in  (4.1)  with  E  [c2]  —  oo.  Then,  with  {£r,  r 
1,2, . . .}  satisfying 


lira 

r— >oo  CL 
Sr 


r  h(Cr)  =  hrn  -~E 

VS  '  r->oo£2 


1  [<7  <  Cr]  O2  —  K 


(4-4) 


for  some  positive  constant  K,  it  holds  that 


{Sr(t),  t>  0} 


B(t),  t>  0}. 


Finally,  we  turn  to  the  case  of  bona  fide  long-range  dependence,  i.e.,  1  <  a  <  2. 
Recall  that  a  Levy  process  is  an  IR— valued  process  with  stationary  independent 
increments.  We  let  {La(t),  t  0}  denote  a  standard  ot— stable,  spectrally  positive 
Levy  motion,  i.e.,  a  Levy  process  such  that  La( 0)  =  0  a.s.  and  for  t  >  0  the  rv 
La{t)  is  a  stable  rv  Sa(f1/a,l,0)  [28]  characterized  by 

E[exp(t0La(i))]  =  exp  (l  -  isgn(0)  tan(™)))  ,  0  6lR.  (4.5) 

Theorem  4.3  (Long-range  dependence)  If  1  <  a  <  2  in  (4.1),  then  with  {Cr,r  = 
1,2,...}  satisfying 


lim  T-h(Cr)  =  lim  rP  [a  >  (r\  =  K  (4.6) 

r— >00  fa  r— foo  1 

Sr 

for  some  positive  constant  K ,  it  holds  that 

t>  0}  =^r  {C(K)  La(t),  t  >  0}  (4.7) 


where 


C(K)  = 


f  cKF(2  —  a) 

l,(a-l)E[<7] 


COs(7T 


2  —  a 


>) 


1/q 


(4.8) 


We  close  with  a  characterization  of  the  scaling  sequences  encountered  in  Theorems 
4.2  and  4.3;  its  proof  is  available  in  Proposition  10.3  of  Section  10. 


Proposition  4.4  The  scaling  sequence  {Cr,  r  =  1,2, ...}  of  Theorems  4.2  and  4.3 
is  1  / a-regularly  varying,  1  <  a  <  2,  i.e.,  of  the  form  (r  =  rl,ah(r)  for  some  slowly 
varying  function  h  :  11+  -»  1R+ . 


5  Consequences  and  comments 

Several  interesting  inferences  follow  from  the  heavy  traffic  results  obtained  so  far. 


11 


5.1  Queue  size 

We  start  with  the  heavy  traffic  behavior  of  the  normalized  queue  length.  Whenever 
the  convergence  (3.10)  holds,  we  can  immediately  conclude  from  (3.9)  and  from  the 
continuity  of  the  reflection  mapping  (via  the  continuous  mapping  theorem  [4,  Thm. 
5.1,  p.  30])  that 

{Qr(t),  t>0}=>r  {Q(t),  t  >  0}  (5.1) 

with 

Q(t)  =  S(t)  -7t-  inf  t(S(x)  -  jx) ,  t  >  0.  (5.2) 

The  form  of  the  limit  derives  from  (3.9)  and  the  fact  that  ^limp^f)  =  —7 1  under 
Assumption  (A). 

This  observation  can  now  be  used  to  provide  a  characterization  of  the  steady- 
state  buffer  content  in  heavy  traffic  under  the  assumptions  of  Theorems  4.1-4.3. 

In  the  short-range  dependent  case,  Theorem  4.1  combines  with  a  classical  result 
on  the  supremum  functional  of  Brownian  motion  [12,  p.15]  to  yield  the  following. 

Theorem  5.1  Under  the  assumptions  of  Theorem  4.1,  the  resulting  stationary 
heavy-traffic  buffer  content  is  exponentially  distributed  with 

Um  P  [0(f)  >  0  =  «P  (-  ")  •  "SO- 

Theorem  4.2  leads  via  (5.1)-(5.2)  to  a  similar  result. 

Theorem  5.2  Under  the  assumptions  of  Theorem  4.2,  the  resulting  stationary 
heavy-traffic  buffer  content  is  exponentially  distributed,  with 

fim  P  [Q{t)  >  b]  =  exp  ^  -  bj  ,  b  >  0. 

Finally,  in  the  stable  case  we  combine  Theorem  4.3  with  the  developments  in  [5, 
Theorem  12a]  to  get  the  following  fact. 

Theorem  5.3  Under  the  assumptions  of  Theorem  4.3,  the  distribution  of  the  re¬ 
sulting  stationary  heavy-traffic  buffer  content  has  a  Pareto  tail  given  by 

lim  P  [Q(t)  >  b]  ~  — - WyFTTT  bl~°  ( b  °°)-  t5’3) 

t->o o  L  7(a-l)E[aJ 
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5.2  On  selecting  the  heavy  traffic  scaling 

As  the  appropriate  scaling  sequence  {(r,  r  =  1,2, . . .}  is  revealing  of  the  nature  of 
the  limiting  heavy  traffic  process  {S{t),  t  >  0},  we  briefly  discuss  here  its  selection. 

It  is  known  [16]  that  convergence  of  a  normalized  partial  sum  process,  such  as 
{Sr(t),  t  >  0},  can  only  be  to  a  self-similar  process,  and  that  the  correspond¬ 
ing  Hurst  parameter  H  may  be  determined  through  the  regularly  varying  scaling 
{0,  r  =  1,2,...}  by 

lim  =  xH ,  x  >  0. 
r->  oo 

On  the  other  hand,  under  (4.1)  with  1  <  a  <  2,  the  M\GI\oo  busy  server  pro¬ 
cess  process  already  possesses  the  so-called  second  order  asymptotic  self-similarity 
property,  with  parameter  (3  -  a)/ 2  [18],  i.e.,  by  aggregating  the  original  process 
{bn,  n  =  0, 1, . . .}  in  blocks  of  size  m  and  dividing  by  the  block  size,  we  can  obtain 
in  the  limit  as  m  goes  to  infinity  the  same  correlation  function  as  that  of  a  fractional 
Brownian  motion.  Because  of  this  fact  it  would  be  tempting  to  think  that  the  ap¬ 
propriate  scaling  ensuring  (3.10)  might  be  the  one  that  balances  the  rate  of  growth 
of  the  partial  sums  variance,  so  that  convergence  occurs  to  a  limiting  process  with 
finite  variance.  By  standard  calculations,  we  find  the  partial  sums  variance  to  be 


for  all  r  =  1, 2, . . ..  As  shown  in  [32],  when  the  tail  of  a  satisfies  (4.1)  with  1  <  a  <  2, 
the  candidate  scaling  {(r,  r  =  1,  2, . . .}  given  by 

C2  =  r^2  p  [?  >  k] ,  r  =  1,2,...  (5.5) 

k= l 


var  [S,r(f)]  = 


ArE  [ct] 


M 

'rt\ +2Y, 


indeed  results  in  a  finite  limiting  variance,  i.e.,  lim  var  [5r(i)]  exists  and  is  finite  for 

r  — >oo 

all  t  >  0.  In  addition,  invoking  (9.24)  we  see  that  this  scaling  has  the  asymptotic 


r2  _ i _  3 

~  (2  —  a)(a  —  1)E  [a] 


*h(r)  (r  -»  oo). 


As  such  a  scaling  is  clearly  (3  -  a) /2-regularly  varying,  it  suggests  possible  con¬ 
vergence  to  a  fractional  Brownian  motion  with  Hurst  parameter  (3  —  a)/ 2.  For 
a  single  M\GI\oo  process  however  this  is  not  true  as  such  convergence  to  a  frac¬ 
tional  Brownian  motion  does  not  take  place.  In  fact,  the  candidate  scaling  (5.5), 


which  balances  the  growth  of  the  variance,  is  too  strong  and  yields  convergence  to 
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a  degenerate  limit  -  the  identically  zero  process.  Theorem  4.3,  in  conjuction  with 
Proposition  4.4,  clearly  shows  that  the  correct  scaling  does  not  contain  any  r^3  a^2 
factor;  but  instead  contains  the  weaker  r1/0  factor  associated  with  the  stable  law  to 
which  the  service  rv  a  is  attracted.  As  a  result,  the  limiting  heavy  traffic  process 
turns  out  to  be  not  a  fractional  Brownian  motion  but  an  a-stable  1 /a  self  similar 
Levy  motion,  the  stable  analog  of  standard  Brownian  motion,  which  has  infinite  vari¬ 
ance  and  independent  increments.  In  heavy  traffic,  the  corresponding  queue  length 
asymptotics  are  not  Weibullian,  but  hyperbolic  in  nature,  with  power  1  -  a.  Thus, 
the  M\GI\oo  processes  demonstrate  that,  within  long-range  dependence,  fractional 
Brownian  motion  does  not  assume  the  ubiquitous  role  that  its  short-range  depen¬ 
dent  counterpart,  standard  Brownian  motion,  plays  in  the  short-range  dependence 
setup,  and  that  modeling  possibilities  attracted  to  non-Gaussian  limits  are  not  so 
hard  to  find.  Clearly,  the  extent  to  which  such  non-Gaussian  processes  can  serve  as 
useful  traffic  models  deserves  some  further  consideration. 

6  Outline  of  proof  and  preliminary  results 

In  this  section  we  organize  the  proof  of  Theorems  4.1  4.3  into  a  series  of  steps  which 
we  formalize  as  Propositions;  their  proofs  are  given  in  Section  8. 

Look  at  the  rth  queueing  system,  r  =  1,2,...,  and  fix  n  =  0, 1,  —  We  note  the 
decomposition  brn  =  b{n)r  +  b{n)r  where  the  rvs  ^°)r  and  bn)r  describe  the  contribu¬ 
tions  to  the  number  of  customers  in  the  system  at  the  beginning  of  slot  [n,  n  +  1) 
from  those  initially  present  (at  n  =  0)  and  from  the  new  arrivals,  respectively.  It  is 
easy  to  see  that 

r  ( 3r 

b£)r  =  E  1  >  n]  and  hn)r  =  E  E  1  Kj  >n-k\. 

j= l  fc=lj=l 

It  was  shown  in  [22]  that 

40)r  =  E  =  X]  min(n>  °j  -  !)  (6d) 

j= i  j= i 

alld  n  n  ft 

4a)r  =  E  4a)r  =  E  E  n  -  k  +  1).  (6.2) 

k= 1  fc=lj=l 
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We  introduce  the  rescaled  versions 


s<o)rw = t  <-gsr 


E 


„(°)r 


D 


and  SMr(t)  3  i  (sj“)|r 


E 


Sa)r 


])• 


t  >  0 


so  that 

Sr(t)  =  SWT(t)  +  S^a)r(t),  t  >  0.  (6.3) 

Also,  for  each  T  >  0,  the  identically  zero  mapping  on  [0,  T]  is  the  element  of  D[ 0,  T] 
denoted  by  9t ,  he.,  6t  '■  [0,T]  — >  1R  with  9x{t)  =  0,  0  <  t  <  T. 

We  first  show  that  the  initial  condition  plays  no  role  in  the  heavy  traffic  limit, 
as  should  be  expected.  This  reduction  step,  as  well  as  others  taken  in  this  section, 
is  accomplished  under  the  following  sufficient  condition. 


Condition  (B)  The  scaling  sequence  {Cr,  r  =  1,2,...}  satisfies 


lim  P  [B  >  j] 
™  Cr  ^  1 


0. 


Condition  (B)  holds  under  each  set  of  assumptions  of  Theorems  4. 1-4. 3;  this  is 
shown  in  Proposition  7. 1  of  Section  7. 


Proposition  6.1  Under  Condition  (B),  for  each  T  >  0  it  holds  that 

{S^r(t),  0  <  t  <  T)  =^r  9t  in  D[0,T}. 

Thus,  in  order  to  get  (3.10)  it  suffices  to  consider  the  limiting  behavior  of  the 
rescaled  process  {5^a)r(t),  t  >  0}.  To  that  end,  for  each  r  =  1,2,...,  we  introduce 
the  sequence  {w^,  n  =  0, 1, . . .}  given  by 

f y 

wro  =  0,  wh  =  Y.J2  n=l,2,...  (6.4) 

k=\ j=i 

which  can  be  interpreted  as  the  contribution  to  the  workload  due  to  arrivals  only. 
The  corresponding  rescaled  process  {Wr(t),  t  >  0}  is  given  by 

wr(t)  =  4  («.[.„ -E[«fJ),  f>0.  (6.5) 

Sr 

The  main  idea  driving  the  discussion  is  that  in  as  much  as  heavy  traffic  is  concerned, 
the  process  {Wr(t),  t  >  0}  acts  as  a  surrogate  for  {S^r(t),  t  >  0}.  This  is  made 
precise  through  the  following  “convergence  together”  result. 


15 


Proposition  6.2  Under  Condition  (B),  for  each  T  >  0  it  holds  that 

{Wr(t)  -  S{a)r{t),  0  <  t  <  T}  =^r  eT  in  D[0,T\. 

Combining  Propositions  6.1  and  6.2,  we  immediately  get  the  following  conclusion 
from  (6.3). 

Corollary  6.3  Under  Condition  (B),  for  each  T  >  0  it  holds  that 

{ Wr(t )  -  Sr(t),  0  <t<T}=^r0T  in  D[0,T], 

so  that  the  process  {Sr{t),  0  <  t  <  T]  is  S-weakly  convergent  if  and  only  if 
{Wr(t),  0  <t<T}is  S-weakly  convergent,  and  convergence  is  to  the  same  limit. 

Thus,  we  need  only  consider  the  convergence  of  the  process  {WT(t),  t>  0},  and 
characterize  the  limiting  process.  In  fact,  a  further  reduction  can  be  achieved  by 
noting  that  in  heavy  traffic  we  can  replace  {(drk,k  —  1,2,...}  by  the  limiting  i.i.d. 
sequence  {/3k,k  =  1,2,...},  where  the  generic  rv  /?  is  a  Poisson  rv  with  parameter 
c/E  [<r].  More  precisely,  consider  the  modified  workload  process  {vn,  n  =  0, 1, . . .} 
given  by 

n  t>k 

vo  —  0;  n  =  1,2,...  (6.6) 

k=lj=l 

under  the  assumption  that  the  rvs  {/3k,  k  =  1,2,...}  are  independent  of  the  service 
time  rvs  {onJ,  n,j  =  1,2,...}.  For  each  r  =  1,2,...,  the  corresponding  rescaled 
process  { Vr(t ),  t  >  0}  is  defined  by 

Pr(f)  =  i(u[rt]-E[u[rt]]),  t>  0.  (6.7) 

Proposition  6.4  Under  Assumption  (A),  the  process  {Wr(t),  0  <  t,  <  T}  is  S- 
weakly  convergent  if  and  only  if  {Vr(t),  0  <  t  <  T}  is  S-weakly  convergent,  and 
convergence  is  to  the  same  limit. 

Corollary  6.3  and  Proposition  6.4  together  lead  to  the  following  conclusion: 

Corollary  6.5  Under  Assumption  (A)  and  Condition  (B),  the  process  {Sr(t),  0  < 
t  <  T\  is  S-weakly  convergent  if  and  only  if  {Vr(t),  0  <  t  <  T}  is  S-weakly 
convergent,  and  convergence  is  to  the  same  limit. 
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7  Proofs  of  Theorems  4. 1-4. 3 

First,  the  big  picture:  Corollary  6.5  and  Proposition  7.1  (given  below)  imply  that 
in  proving  Theorems  4. 1-4.3  we  need  only  investigate  the  convergence  of  the  mod¬ 
ified  workload  process  (6.7).  This  is  a  much  easier  task  as  we  now  deal  with  the 
(normalized)  partial  sums  process  associated  with  a  single  sequence  of  i.i.d.  rvs, 
of  finite  mean  but  possibly  infinite  variance,  an  extensively  studied  situation  where 
the  (functional  form  of  the)  classical  Central  Limit  Theorem  and  its  generalization 
to  i.i.d.  summands  with  infinite  variance,  are  expected  to  yield  the  requested  con¬ 
vergence.  In  fact,  as  we  shall  see  shortly,  the  convergence  of  the  finite  dimensional 
distributions  of  {Vr{t),  t  >  0}  turns  out  to  be  an  easy  by-product  of  classical 
results  concerning  stable  distributions  and  their  domains  of  attraction  [11,  pp.  574- 
581],  Finally,  the  desired  S-weak  convergence  of  the  process  {Vr(t),  t  >  0},  thus 
of  {Sr(t),  t  >  0},  will  be  validated  through  the  functional  convergence  results  due 
to  Skorokhod  [29].  This  approach  clearly  explains  the  form  of  the  results  obtained 
in  this  paper,  providing  insights  as  to  when  the  process  {Vr(t),  t  >  0}  is  expected 
to  converge,  and  to  which  limit. 

And  now,  on  with  the  details:  In  Section  9  we  give  a  proof  that  the  technical 
Condition  (B)  required  to  establish  the  “convergence  together”  argument,  indeed 
holds  under  the  assumptions  of  Theorems  4. 1-4.3. 

Proposition  7.1  Condition  (B)  holds  true  for  each  of  the  scaling  sequences  {Cr,  r  = 
1,  2, . . .}  in  Theorems  4. 1-4.3. 

Next,  we  consider  the  generic  compound  rv  Y  given  by 

¥=•£0,  (7.1) 

3= 1 

where  the  rv  (5  is  a  Poisson  rv  with  parameter  c/E  [cr]  and  independent  of  the  i.i.d. 
rvs  {cry ,  j  =  1,2, . . .}  which  are  distributed  according  to  cr.  Fixing  t  >  0,  we  note 
that 

,  M 

vr(t )  =st  -  £  (r*  -  E  [Yk]) ,  r  —  1,2, . . .  (7.2) 

^  k= 1 

where  the  i.i.d.  rvs  {Tfc,  k  =  1,2,...}  are  distributed  according  to  the  generic  rv  Y . 

For  easy  reference,  we  restate  some  useful  facts  concerning  stable  distributions 
and  their  domains  of  attraction;  the  reader  is  refered  to  [11,  pp.  574-581]  for 
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additional  material:  Let  L  be  a  rv  with  distribution  not  concentrated  at  one  point, 
and  let  {Xr,  r  =  1, 2, . . .}  be  a  sequence  of  i.i.d.  rvs,  with  generic  rv  X.  We  say 
that  X  belongs  to  the  domain  of  attraction  of  the  rv  L  if  there  exist  normalizing 
constants  Cr  >  0  and  cr,  r  =  1,2,...,  such  that 


X\  +  .  .  .  +  Xr  —  TCr 
<; 


L. 


(7.3) 


By  Theorem  1  of  [11,  p.  576]  only  stable  rvs  possess  a  domain  of  attraction.  By 
Theorem  2  in  [11,  p.  577],  in  order  for  X  to  belong  to  the  domain  of  attraction  of 
a  stable  law  with  exponent  a,  0  <  a  <  2,  it  is  necessary  that  its  truncated  second 
moment  be  regularly  varying  with  exponent  2  —  a,  i.e., 


E  1  [X  <  r]  X2]  ~  r2  ag(r)  (r  — >  oo), 


(7.4) 


for  some  slowly  varying  function  g  :  IR+  — >  1R+.  The  associated  scaling  sequence 
{Cr>  r  =  1,2,...}  in  (7.3)  must  satisfy 


lim^E[l[X  <Cr]X2 

r— >  ooQf  L 


=  M 


(7.5) 


for  some  constant  M  >  0  [11,  p.  579].  Moreover,  if  E  [X]  is  finite,  then  by  Theorem 
3(ii)  of  [11,  p.  581]  we  can  take  cr  =  E[X],  r  =  1,2, — 

We  are  now  ready  to  discuss  Theorems  4. 1-4.3  which  are  all  proven  in  the  same 
manner,  although  for  clarity  of  presentation,  we  shall  consider  each  of  them  sepa¬ 
rately.  As  E  [Y]  is  finite  under  the  enforced  assumptions,  we  conclude  from  (7.2)  and 
(7.3)  that  for  each  t  >  0,  the  convergence  question  concerning  (Fr(t),  r  =  1,2,...} 
is  equivalent  to  determining  whether  the  rv  Y  is  attracted  to  a  stable  law,  and  to 
which  one.  In  asserting  this  equivalence  we  rely  on  the  fact  that  the  scaling  sequence 
{Cri  r  ~  1,2,...}  so  selected  is  regularly  varying;  this  turn  out  to  be  the  case  by 
Proposition  4.4,  so  that 

lim  =  t1/a ,  t  >  0.  (7.6) 

r—yoo 

In  each  case,  we  show  that  both  the  necessary  condition  (7.4)  and  the  accompanying 
sufficient  condition  stated  in  [11,  p.  577]  are  satisfied.  This  occurs  simply  because 
the  generic  rv  Y  inherits  the  tail  behavior  of  the  generic  service  time  a  under  each 
set  of  assumptions  of  Theorems  4. 1-4.3. 


A  proof  of  Theorem  4.1.  Since  E  [<r2]  <  oo,  the  variance  of  Y  is  also  finite,  and 
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is  given  by 


(7.7) 


var  [Y]  =  var  [/3\  E  [a]2  +  E  [/?]  var  [er]  =  ° g 

Hence,  the  truncated  second  moment  of  Y  varies  slowly,  i.e. ,  (7.4)  holds  with  a  — 
2  and  as  Y  is  never  degenerate  at  one  point,  it  follows  from  Corollary  1  to  Theorem 
2  in  [11,  p.  578]  that  Y  is  attracted  to  the  normal  distribution.  Obviously,  the 
scaling  (r  =  y/r,  r  —  1,2, .. .  satisfies  (7.5),  with  M  =  cE  [a2]  /E  [a].  In  fact,  by  a 
well-known  result  of  Donsker  [4,  Thm.  16.1,  p.  137],  selecting  £r  =  y/r ,  r  —  1,2,... 
ensures  that  the  process  {Vr(t),  t  >  0}  is  5- weakly  convergent  to  a  Brownian 
motion,  with 

{Vr(t),  t  >  0}  =^r  {v/M  B(t),  t  >  0}.  (7.8) 

Combining  (7.8)  with  Proposition  7.1  and  Corollary  6.5  immediately  concludes  the 
proof.  ■ 

Under  the  assumptions  of  Theorems  4.2  and  4.3,  E  [a2]  is  infinite,  and  the  compound 
Poisson  rv  Y  now  has  infinite  variance.  Also,  if  a  satisfies  the  tail  condition  (4.1), 
so  does  Y  with 


P  [Y  >  r)  =  P 


a1  >  r 


3= i 


E[/0]  r  ah(r)  (r  — »■  oo). 


(7.9) 


The  asymptotic  equality  in  (7.9)  is  stated  as  an  exercise  in  [11,  Ex.  31,  p.  288], 
where  the  reader  will  find  hints  for  its  proof.  Next,  we  check  that  the  truncated 
second  moment  of  Y  is  given  by 


E  [l  [Y  <  r]  Y2 


r  ) — 1 

=  2  ^2  \Y  >  n]  —  r(r  +  2)P  [Y  >  r]  +  y]  P  [Y  >  n ] 

n= 1  n= 0 


for  each  r  =  1,2,...,  and  using  (7.9)  in  this  last  expression  we  find  that 


E  [l  [Y  <  r]  y2]  ~  E  [p]  I  2  J2  nl~ah{n ) 

\  n=l 


(r  — >  oo) 


(7.10) 


because  E  [y]  is  finite  and  E  [y2]  is  infinite.  We  close  these  preliminary  remarks 
by  noting  that  the  truncated  second  moments  of  a  and  Y  are  obviously  related  to 
each  other  by 

E  [l  [y  <  r]  y2]  ~  E  [0\  E  [l  [a  <  r]  a2]  (r oo).  (7.11) 
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A  proof  of  Theorem  4.2.  Inserting  a  =  2  in  (7.10)  and  using  the  definition  (4.2) 
(with  a  —  2),  we  get 

E  [l  [Y  <  r]  y2]  ~  E  [fi\  (h(r)-h(r))  (r  -+  oo). 


By  Proposition  9.2(ii),  l2  :  1R+  ->  K+  is  slowly  varying  with 


lim 

T— >00 


HH 

Hr) 


=  0, 


so  that 

E  [l  [Y  <  r]  y2j  ~  E  \/3\  h(r)  (r -»  oo).  (7-12) 

This  time,  by  Corollary  1  in  [11,  p.  578],  the  slow  variation  of  the  truncated  second 
moment  is  a  necessary  and  sufficient  condition  for  Y  to  be  attracted  to  the  nor¬ 
mal  distribution,  with  normalizing  coefficients  selected  according  to  (7.5)  (despite 
the  fact  that  the  variance  of  Y  is  now  infinite).  Since  the  marginals  of  the  process 
{Vr{t),  t  >  0},  which  has  stationary,  independent  increments,  converge  to  a  Gaus¬ 
sian  distribution,  it  follows  by  [29,  Theorem  2.7]  without  any  additional  conditions 
that  (7.8)  takes  place.  Because  of  (7.12),  selecting  the  scaling  {Cr,  r  =  1,2,...} 
according  to  (7.5),  with  M  =  cK/E  [cr],  is  equivalent  to  (4.4).  Combining  (7.8)  with 
Proposition  7.1  and  Corollary  6.5  completes  the  proof.  ■ 


A  proof  of  Theorem  4.3.  When  1  <  ct  <  2  in  (4.1)  the  rvs  cr  and  Y  have 
infinite  second  moment,  and  Proposition  9.2 (i)  implies 

1 


lim 


7’— >007-2  ah(r) 

v  ’  n— 1 


Y, nl  aHH  = 


(2  —  a) 


(7.13) 


Using  this  asymptotic  in  (7.10)  we  get 

E  |l  [Y  <  r]Y2]  ~  Q-~  E  [f3\  r2~ah{r )  (r oo).  (7.14) 

Invoking  Corollary  2  of  [11,  p.  578],  we  see  that  (7.14)  and  the  tail  condition  (7.9) 
are  sufficient  to  ensure  that  Y  belongs  to  the  domain  of  attraction  of  a  non-normal 
stable  distribution  with  exponent  1  <  a  <  2.  The  associated  scaling  sequence 
{(r,  r  =  1,2,...},  selected  according  to  (7.5),  yields  convergence  of  the  marginal 
distribution  of  Ur(l),  as  r  goes  to  infinity,  to  that  of  an  a-stable  rv,  i.e., 


lim  E  [exp(i#Ur(l))]  =  E 


exp 


iO 


(  MT  (3  —  a) 

V  a(a  —  1) 


cos(7r 


2  —  a 


') 


l/o 


La(  1) 


0  G  1R . 
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The  exact  value  of  the  constant  given  above  can  be  easily  verified,  by  recalling  the 
expression  (4.5)  for  the  characteristic  function  of  La(l)  and  comparing  it  with  Eq. 
(3.18)  of  [11,  p.  730]  (note  the  unfortunate  error  in  the  ±  sign).  Next,  appealing 
to  [29,  Theorem  2.7]  again,  we  conclude  that  convergence  of  the  marginals  also 
implies  5-weak  convergence  of  the  process  \Vr(t),  t  >  0},  which  has  stationary, 
independent  increments,  to  an  a-stable  Levy  motion.  More  precisely,  it  holds  that 

{Vr(t),  t  >  0}  =4r  {(MT{3  ~  cos(t T^))1/a  LQ(t),  t>  0}.  (7.15) 

V  a(ce  -  1)  2  / 

Using  (7.14)  in  (7.5)  with  M  =  cKa/( 2  —  a)E  [a]  we  obtain  the  scaling  sequence 
{£ r ,  r  =  1,2,...}  given  in  (4.6).  Finally,  combining  (7.15)  with  Proposition  7.1  and 
Corollary  6.5  shows  that  (4.7)  holds  true.  ■ 


8  Proofs  of  Propositions  6.1,  6.2  and  6.4 

We  start  by  remarking  that  if  the  sequence  {£r,  r  =  1,2, . . .}  is  regularly  varying 
(as  stated  in  Proposition  4.4),  then  Condition  (B)  also  implies 

i  M 

Km  -£P[a>j]=0,  t>  0.  (8.1) 

c,r  j= i 

All  three  proofs  given  in  this  section  follow  the  same  pattern,  and  are  based  on 
the  following  simple  idea:  Consider  a  sequence  of  IR-valued  processes  {Xr(t),  t  >  0}, 
r  —  1,2,...,  with  sample  paths  in  D[ 0, oo).  Fix  T  >  0.  According  to  Theorem  4.1 
of  [4,  p.  25],  in  order  to  prove  the  U  weak  convergence 

{Xr{t),  0  <t  <T}  6T  in  Z?[0,T],  (8.2) 

it  suffices  to  establish  the  convergence  in  probability 

sup  |Xr(t)|  -A r  0.  (8.3) 

0  <t<T 

A  proof  of  Proposition  6.1.  Fix  r  =  1,2, . . .,  and  note  from  (6.1)  that 

sup  |S'(0)r(f)|  <  min(T?  _  lrT ])  +  Ar~f  —  E  [min(fi  -  1,  [rT])] . 

0 <t<T  Sr  Sr 
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Hence,  for  every  e  >  0,  it  is  plain  that 


P  sup  |5^r(f)|  >  £ 

0  <t<T 

<  P  -3- ^min(CTj  -  1,  [rT])  +  E[min(a  -  1,  [rT])]  >  e 

>r  i= i  Cr 

<  2ArE  [<T]-  E  [min(a  -  1,  [rT])]  (8-4) 

-Cr 

where  the  last  step  follows  by  Chebyshev’s  inequality.  It  is  also  the  case  that 

[rT]-l 

E  [min(<T  -  1,  [rT])]  =  ^  P  [min(cx  -  1,  [rT])  >  n] 

n=0 
[rT\- 1 

=  Y  p  - 1  >  nl 

n— 0 
[rT] 

=  Y p  &  >  n]  • 

71  =  1 

Appealing  to  Condition  (B)  and  (8.1),  we  get 

lim  E  min(<7  —  1,  [rT])  =0  (8.5) 

r->  oo  [Cr 

and  the  conclusion 

lim  P  sup  |5^°^r(t)|  >  £  =0 
r->o°  [o<t<T 

immediately  obtains  from  (3.2)  upon  letting  r  go  to  infinity  in  (8.4).  ■ 


A  proof  of  Proposition  6.2.  Fix  r  =  1,2, . . ..  From  (6.2)  and  (6.4)  we  note  that 

fiT 

<  -  4a)r  =  Y  Y(akj  -  (n  -  k  +  1))+ 

k=lj=l 

or 

n  Pj-k+1 

=  Y  Y  K-ik+ij  -  k)+ 

k=i  j= i 
n  @k 

=st  YY(ak<j~k)+'  n==1’2’"-  (8.6) 

k= 1 j= 1 


22 


where  in  the  last  step  we  used  the  mutual  independence  of  the  families  of  i.i.d.  rvs 
{/ 3rk ,  k  =  1,2,...}  and  {akj,  k,j  =  1,2,...}  are  mutually  independent.  It  is  now 
straightforward  to  check  that 

i  [rTl  Fh  A  [rT] 

sup  \Wr(t)-S^r(t)\  <st-EEKn-fc)+  +  7Z1EE[((J-fc)+]-  (8-7) 

0 <t<T  W=m=i  k= 1 

By  Chebyshev’s  inequality,  for  every  e  >  0  we  obtain 


1  2A  [rT] 

sup  \Wr(t)-S^r(t)\>e  <  — ^  E  [(a  -  k)+] 

><t<T  J  £l>r  k=1 


£  p  [<?  >  k] , 

esr  t=i 


and  the  desired  convergence 

sup  |W'r(i)-S(°)r(i)|  -Ar0  (8.8) 

0  <t<T 

follows  upon  letting  r  go  to  infinity  in  the  upper  bound  (8.8),  and  making  use  of 
(3.2),  Condition  (B)  and  (8.1).  ■ 


The  proof  of  Proposition  6.4  requires  estimates  that  derive  from  various  mar¬ 
tingales  inequalities;  we  now  state  them  in  Lemmas  8.1  and  8.2  for  easy  reference: 
Consider  a  collection  of  integrable  rvs  {Xu  i  =  1, . . . ,  n}  adapted  with  respect  to  the 
filtration  {Ti,  i  =  1, . . .  ,n},  i.e.,  for  each  i  =  1, . . . ,  n,  the  rv  Xi  is  .^-measurable. 
We  also  write 

Sl  =  Xl  +  ...  +  Xl,  i  =  1, . . .  ,n.  (8.9) 

Lemma  8.1  (Maximal  inequality  [13])  Assume  {(5,,^)),  i  =  l,...,n}  to  form  a 
martingale.  Then,  for  each  p  >  1,  it  holds  that 

P  max  | S'* |  >  A  <  A~PE  [|Sn|p] ,  A  >  0. 

Lemma  8.2  (von  Bahr  -  Esseen  inequality  [2])  Assume  {(Xi,Xi),  i  =  1, . . . ,  n} 
to  form  a  martingale  difference.  If  E  [|2C|P]  <  oo  for  all  i  =  1, . . . ,  n,  then 

E[|Snn<2^E[|XlH,  1<P<2. 

i= 1 
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In  what  follows  Lemmas  8.1  and  8.2  are  applied  to  the  case  when  the  rvs  {X,,  i  = 
1, . . . ,  n}  are  zero-mean  i.i.d.  rvs. 


A  proof  of  Proposition  6.4.  Recall  that  the  rvs  {/3k,  k  =  1,2,...}  are  i.i.d. 
Poisson  rvs  with  parameter  c/E[cr],  which  are  independent  of  the  sequence  of  i.i.d. 
service  rvs  {crkj,  k,j  =  1,2, . . .}. 

Fix  r  =  1,2,....  On  the  same  probability  triple  (Q,X,  P)  where  the  previ¬ 
ously  mentioned  rvs  are  defined,  we  introduce  a  family  of  i.i.d.  {0,l}-valued  rvs 

{UrkJ,  k,j  =  1,2,...},  i.e., 


P  [Ur  =  1]  =  =  I  _  p  [Ur  =  o] 


where  Ur  denotes  the  generic  rv  for  this  i.i.d.  sequence.  The  rvs  {Uk  j,  k,j  = 
1,  2, . . .}  are  assumed  independent  of  the  collections  of  rvs  mentioned  so  far.  Next, 
we  define  the  rvs  {/3rk,  k  =  1, 2, . . .}  by 

j= i 


We  also  define  the  workload  process  { w Tn,  n  =  0, 1, . . .}  corresponding  to  {/?£,  k  = 
1,2,  ...}by 

n  &r 

wT0  =  0,  n=  1,2,... 

k= li=l 

and  its  rescaled  version  {Wr(t),  t  >  0}  by 


W'(t)  =  ±-(w'[rt]-E  [»frt|]),  (>0. 


Under  the  enforced  independence  assumptions,  it  is  easy  to  check  that  {/3},  k  = 
1,2,...}  =st  {/3rk,  k  =  1,2,...},  and  that 


{Wr(t),  0  <  t  <  T}  =st  {Wr(t),  0  <  t  <  T}. 


Moreover,  these  rvs  are  all  defined  on  the  same  probability  triple  as  the  rescaled 
process  {Vr(t),  t  >  0}.  Thus,  the  result  will  be  established  if  it  holds  that 

{Vr{t)  -  Wr(t),  0  <  t  <  T}  0T  in  D[0,T], 
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or  equivalently,  if  we  can  show  that 


sup  \Vr(t)  -  Wr{t)\  -Ar  0.  (8.10) 

0  <t<T 

To  that  end,  for  each  r  =  1, 2, . . we  note  from  the  definitions  that 

n  0k 

^  =  n  =  l,2,... 

j=i 


so  that 


n  0k 

n  ~  Ukj)°kj,  n=  1,2,... 

fc=l j=l 


r  ~r 

vl  -  W 


The  rvs  {Zrk,  k  =  1,  2, . . .}  defined  by 


0k 

Zl  =  £(1  -  Ul3)ak,v  k  =  1,2, 
l=i 


(8.11) 


(8.12) 


are  i.i.d.,  and  we  denote  by  Zr  the  corresponding  generic  rv  associated  with  this 
collection  of  rvs.  It  is  plain  from  (8.11)  and  (8.12)  that 


sup  \Vr{t)  -  Wr(t)\  =  ~  sup 

0 <t<T  Sr  l<n<[rT] 


£(ZJ-E[ZJ]) 


fc=l 


(8.13) 


Fix  e  >  0.  Invoking  the  maximal  inequality  for  martingale  sequences  [Lemma  8.1], 
we  get 


P 


sup  I Vr(t)  -  Wr{t) |  >  £ 
0  <t<T 


1 

■  RTjp 


E 


[rT] 

E^-E[Z[]) 

fc= l 


(8.14) 


with  p  selected  such  that  1  <  p  <  a  <  2.  This  selection  of  p  ensures  E  [ap]  <  oo 
both  under  short-range  dependence  and  under  the  assumption  of  regularly  varying 
tail  (4.1).  The  von  Bahr  -  Esseen  inequality  [Lemma  8.2]  for  martingale  differences 
can  now  be  applied  to  the  right  handside  of  (8.14)  to  yield 


1 

«r)P 


E 


£(Z£-E[Z[]) 

k=l 


<^E{\z'-E{znn- 


By  the  convexity  of  x  — »  xp  (p  >  1)  on  IR+,  we  find 


E  [\Zr  -  E  [Zr]  |p]  <  2P_1(E  [|Zr|p]  +  E  [\Zr\]p)  <  2PV  [|Zr|p] 


(8.15) 


(8.16) 
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with  the  last  step  validated  by  Jensen’s  inequality.  Next,  using  the  definition  of  Zr , 
we  obtain  by  the  same  convexity  argument  that 


E  [|Zr|p|/3]  <  (3p~l  E 


P 

E(i  -  u?P(jPj  w 


U= 1 


=  /3pE[op]E[(l  -  Ur)p] 


a.s.  (8.17) 


under  the  enforced  independence  assumptions  (and  with  an  obvious  notation).  In¬ 
jecting  the  bounds  (8.16)  and  (8.17)  into  (8.15),  we  conclude  from  (8.14)  that 


P 


sup  \Vr(t)  -  Wr{t)\  >  e 

0  <t<T 


<  J^i  [ap]  E  \pp]  ^E  [(1  -  Ur)p } 


(8.18) 


As  the  heavy  traffic  Assumption  (A)  implies 

lim  — E  [(1  —  Ur)p]  =  — ,  (8.19) 

r— Kx>£r  c 

the  desired  conclusion  (8.10)  now  follows  by  letting  r  go  to  infinity  in  (8.18),  using 
(8.19)  and  noting  that  lim  1/(|?-1  =  0  for  p  >  1.  ■ 


9  A  proof  of  Proposition  7.1 

In  the  proof  of  Proposition  7.1  and  elsewhere,  we  make  use  of  the  following  fact. 


Lemma  9.1 

For  any  slowly  varying  function 

u  :  IR+  — >  1R+,  it  holds  that 

lim  xpu(x)  =  oo, 

p>  o, 

(9.1) 

x — >oo 

while 

lim  xpu(x)  =  0, 

X — KXD 

p  <  0. 

(9.2) 

Proof.  By  the  Representation  Theorem  for  slowly  varying  functions  [6,  Theorem 
1.3.1,  p.  12],  we  can  write 

u(x)  ~  cexp  (/  dt'j  ( x  — >  oo)  (9.3) 

with  constants  A  >  0  and  c  >  0,  and  Borel  mapping  e  :  1R+  -»  1R  such  that 
lim  e(t)  =  0.  Thus, 

t— KX> 

xpu(x)  ~  cA~p  exp  (^J  ^  ■-  dt^  (x  -¥  oo).  (9-4) 
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(9.5) 


For  every  5  >  0  there  exists  tg  >  A  such  that  ]£(t)|  <  <5  for  t  >  t, 5,  whence 


-5+p<m±i<s±p 

t  ~  t  -  t 


so  that 


K  +  (—<5  +  p)  In  ^ ^  <  J  - - — —  dt  <  K  +  (5  +  p)  In  )  ,  x  >  tg  (9-6) 


with 


e(t)  +  p 

fts  e(t)  +  p 


K  = 


-r 

.m 


dt. 


I A  t 

The  conclusion  (9.1)  (resp.  (9.2))  immediately  follows  from  these  inequalities  when 
selecting  6  >  0  such  that  5  <  p  (resp.  6  <  -p),  such  a  selection  is  always  possible 
when  p  >  0  (resp.  p  <  0).  ■ 


The  limit  (9.1)  is  useful  in  the  proof  of  the  following  discrete  analogue  to  the 
direct  half  of  Karamata’s  Theorem  [6,  p.  26]. 

Proposition  9.2  Let  u  :  IR+  ->  1R+  be  a  slowly  varying  function.  Then  the  fol¬ 
lowing  statements  hold: 

(i)  For  any  p  >  -1,  we  have  the  asymptotics 

T  rp+l 

^2  npu(n)  ~  -  u(r)  (r  00).  (9.7) 


n— 1 


p  + 1 


(ii)  The  mapping  u  :  IR+  — »  fft+  defined  by 


u 


ix  /  \ 

n= 1 


is  a  slowly  varying  function  which  satisfies 

u(x ) 


lim  '  ~  =  0. 

x-+oo  u(x) 


(9.8) 


(9.9) 


Proof,  (i)  Pick  5  in  (0, 1  +  p).  By  Potter’s  bound  [6,  p.  25],  for  every  A  >  0  there 
exists  rg  a  such  that 


u(n )  ,  fn 


,  .  <  A  -  ,  rgjA  <  n  <  r. 

u(r)  \r 


(9.10) 


27 


Because  p  +  1  >  0,  we  readily  see  from  (9.1)  that 
1 


rP+1u(r)  , 
'  '  n=l 


npu(n ) 


^  1  f  n\p  u(n)  .  . 

E  -  -  -ZT  (r  °°) 

A.  ^  ly  \  r  '  “  1 


n=rs<A 


r  \r  J  u(r) 


- 

/  Ur(x)  dx 

Jo 

(9.11) 

where  for  each  r  =  1,2,... 

,  we  have  defined 

the  functions  Ur  : 

1R+  — ►  1R  and 

Tr  :  1R+  — >  1R  by 

Ur(x)  =  1 

.  rg,A' 
x  >  — — 
r 

Tr(xF  U(r™],  X  >  0 
u{r) 

(9.12) 

and 

ry 

x  >  0. 

(9.13) 

For  every  z  in  (0, 1]  and  r  =  1, 2, . . .  we  note  that  x  <  Tr{x)  <  2,  and  the  Uniform 
Convergence  Theorem  for  slowly  varying  functions  [6,  Theorem  1.2.1,  p.  6]  thus 
implies 

_ .  ( —nn  /  \ 

(9.14) 


u(rTr(x)) 

lim  v  ,rv.  =1,  0  <  x  <  1. 

r—too  u(r) 


Combining  (9.12),  (9.13)  and  (9.14)  we  readily  conclude  to  the  pointwise  convergence 

lim  Ur(x)  =  xp,  0<x<l.  (9.15) 

r— >oo 

Simple  bounding  arguments  based  on  (9.10)  and  the  classical  asymptotics 

1  (  n\p~5 


oo 


lim  E  = - 3-T< 

r^°°n£?SA  r  \rj  p-6  +  1 

can  be  used  to  validate  the  use  of  the  dominated  convergence  theorem,  so  that 

lim  f  Ur(x)dx  =  [  xvdx  =  - -. 

J0  Jo  P  +  1 

The  conclusion  (9.7)  follows  readily  from  (9.11). 

(ii)  Pick  e  in  (0, 1).  For  each  x  >  0  we  have 

r_i 

u(n) 


u(x)  >  E 


n={ex]+l 


n 


=  u{x)  E 


ri=[ex]  +  l 


1  n(x|) 
n  u(x) 


1  j  i 

~  u(x)  E  —  (z  ->  oo), 


(9.16) 


n=[ex}+ 1 
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where  the  asymptotic  equality  follows  by  the  Uniform  Convergence  Theorem  [6, 
Theorem  1.2.1,  p.  6].  Noting  that 


lim  V  —  =  —  In  { 

x—>oo  * — '  fl 

7i= [ex]  4- 1 


(9.17) 


where  e  can  be  chosen  arbitrarily  small,  we  obtain  from  (9.16)  that 

r  "(*)  _ 

lim  -  =  oo 

x^cou{X) 


(9.18) 


or  equivalently,  (9.9). 

To  prove  that  u  is  slowly  varying,  pick  y  >  1  and  note  for  every  x  >  0  that 


l(yx)  =  u{x )  +  u(x)  y 


n=[x]+l 


1  u(n) 
n  u(x) 


Applying  the  Uniform  Convergence  Theorem  [6,  Theorem  1.2.1,  p.  6]  once  more  we 
get 

[P  1  u(n)  1  u(x$)  ...  ,  . 

y - y - -z-~u(x)lny  (i-ioo), 

*ri  ,  n  u(x)  ,  n  u(x) 

n=[x]  +  l  v  '  n=[x]+l 


whence 


u(yx)  ,  u(x)  .  ,  .  , 

/  ~  1  +  ^44  lny  ~  1  (x  -4  oo) 
u(:r)  iqrr) 


upon  using  (9.9).  The  case  y  <  1  is  handled  in  a  similar  way,  and  the  slow  variation 


of  u  follows. 


A  proof  of  Proposition  7.1  From  (2.1)  it  always  holds  that 


~  E[cr2]  1 

B  [5-1  =  x:  P  [^  >  "]  =  2^F]  +  2  ■ 

n= 0  1  J 


(9.19) 


We  consider  each  of  the  scalings  {(T,  r  =  1,2, . . .}  associated  with  Theorems  4.1 
4.3,  separately: 

[Theorem  4.1]  Under  short-range  dependence,  we  have  E  [a2]  <  oo,  and  it  is 
immediate  from  (9.19)  that  Condition  (B)  holds  for  the  choice  Cr  =  y/f,  r  =  1, 2, . . . 
(in  fact  for  any  choice  such  that  lim  (r  =  oo).  ■ 

v  J  r-koo 
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We  next  turn  to  Theorems  4.2  and  4.3.  Upon  substituting  (4.1)  (with  1  <  a  <  2) 
into  (2.1),  we  readily  get  from  Proposition  9.2 (i)  that 

1  00  1 

P  [S  >  nl  =  EH  ?  P  [<r  >  i]  ~  (a-l)'EW  "  h{n)  ("  -»  (9'20) 

J  3 — 71 

whence 

^  P  [d  >  n]  -  - - 1  7  t  nl~ahW  (r  °°)  (9-21) 

(a  —  1)E  [ct] 

provided  E  \5\  is  infinite. 

[Theorem  4.2]  When  a  —  2  in  (4.1),  the  condition  E  [cr2]  =  oo  implies  that  E  [a] 
is  infinite  by  (9.19).  Thus, 

E  P[?>n]~^jEv  (r->oo)  (9.22) 

n=l  f  J  n=l 


which,  from  Proposition  9.2(ii)  is  seen  to  be  slowly  varying.  By  Proposition  4.4,  the 
scaling  {(r,  r  —  1,  2, . . .}  is  1/2-regularly  varying,  so  that 


-T  E  p  >  n]  ~ 

^  n=l 


EHr 


h(r) 


E 

n=l 


h(n) 

n 


( r  -4  oo) 


(9.23) 


for  some  slowly  varying  function  h  :  1R+  -4  1R+.  The  ratio  of  slowly  varying  functions 
being  itself  slowly  varying,  we  readily  conclude  from  Lemma  9.1  and  (9.23)  that 
Condition  (B)  holds.  ■ 


[Theorem  4.3]  On  the  range  1  <  a  <  2,  E  [cr2]  is  infinite,  and  so  is  E  [a]  by  virtue  of 
(9.19).  Proposition  9.2(i)  applied  to  the  right  handside  of  the  asymptotic  equivalence 
(9.21)  yields 


Ep 

71=1 


a  >  n  ~ 


r2  ah(r )  (r  -4  oo). 


(2  -  a)(a  -  l)E[cr] 


(9.24) 


By  Proposition  4.4  the  scaling  {^r,  r  =  1,2, . . .}  is  1 /c*-regularly  varying,  so  that 


y-  E  p  [?  >  nl  ~ 

C’P  71  =  1 


1  h(r)  r2  “ 

(2  -  a) {a  -  1)E  [cr]  h{r)  r« 


(r  -4  oo) 


(9.25) 


for  some  slowly  varying  function  h  :  1R+  — >  1R+.  The  ratio  of  slowly  varying  func¬ 
tions  is  itself  slowly  varying,  and  Condition  (B)  is  now  a  direct  consequence  of 
Lemma  9.1  once  we  note  that  2  —  a  —  a-1  <  0  for  a  >  0.  ■ 
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10  Appendix  -  Asymptotic  invertibility  of  regularly  vary¬ 
ing  functions 


The  proof  of  Proposition  4.4  relies  on  the  following  general  fact. 


Proposition  10.1  Consider  a  Borel  measurable  function  (3  :  1R+  ->  H  such  that 
lim  3{x)  =  0.  With  A  >  0,  the  equation 

x— »oo 


y  =  exp  ^ J  - ^  ^  dtj  ,  x  >  A,y  >  0  (10-1) 

has  a  unique  solution  x  =  x(y)  for  all  y  large  enough.  Moreover, 


lim 

y-yoc 


At  y) 

Ay) 


=  7,  7  >  0, 


(10.2) 


or  equivalently  the  mapping  y  ->  x(y)  is  regularly  varying  of  order  1,  i.e.,  x(y)  ~ 
yv(y)  (  y  — >  oc )  for  some  slowly  varying  function  v  :  1R+  -»  ]R+- 


Proof.  Set 


and  pick  e  in  (0, 1). 


B(x)  =  [X  dt ,  x>  A  (10.3) 

Ja  t 

Since  lim  3{x)  =  0,  there  exists  A*  =  A*{e)  >  A  such  that 
£->00 


l_-_e  l-m  1  +  1  t  >  (10.4) 

ttt 

It  is  straightforward  to  see  that  lim^S(a;)  =  oo  and  B(A)  =  0,  and  by  continuity, 
the  range  of  x  ->  B(x)  contains  the  semi-infinite  interval  [l,oo).  We  also  conclude 
from  (10.4)  that  x  -»  B(x)  is  strictly  monotone  increasing  on  the  interval  [A*,  oo), 
and  the  existence  and  uniqueness  of  a  solution  to  (10.1)  follows  whenever  y  >  y* 
with  y*  =  exp (B(A*)).  The  solution  mapping  y  ->  x{y)  is  strictly  increasing  on 

[y*,  oo). 

We  now  turn  to  proving  (10.2).  There  is  nothing  to  prove  when  7  =  1.  With 
7  >  1,  whenever  y  >y*,  we  get 

—  =  exp(£(x(7 y))  -  B{x{y)))  =  exp  (  f  - — dt  \  ,  (10.5) 

y  \Jx(y)  1  j 


and  the  use  of  the  inequalities  (10.4)  yields 


'ArryY\l 

.  Ay) . 


<  7  < 


'ArryY 
.  Ay) . 


l+£ 


(10.6) 
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or  equivalently, 


(10.7) 


_j_  x('yy)  _  _j_ 

'■yl+e  <C  -  <  'y  1-c  , 


x{y) 


Letting  y  go  to  infinity  in  (10.7),  we  get 


i  .  „x{yu)  x('YV)  _i_ 

7  <  lim  inf  — — —  <  lim  sup  — r—  <  7 1~c , 

1/^00  a;(y)  x{y) 


(10.8) 


and  (10.2)  is  obtained  as  we  note  that  e  is  arbitrary  in  (0, 1).  The  case  7  <  1  is 
handled  in  a  similar  way;  details  are  omitted  in  the  interest  of  brevity.  ■ 


Lemma  10.2  Consider  slowly  varying  functions  u,  w  :  1R+  — >  M+,  such  that  u(x)  ~ 
w(x )  ( x  ->  00),  and  let  a  >  1.  For  any  sequences  {(r,  r  =  1,2,...}  and  {rjr,  r  = 
1,2,...}  with  lim  Cr  =  lim  (r  =  00  such  that 

r— >00  r— >•  00 

lim  r(~au( (r)  =  lim  rrj~aw(r]r)  =  K  (10.9) 

r— >00  r— kx) 

for  some  finite  constant  K  >  0,  it  holds  that  Cr  ~  rjr  (r  00). 


Proof.  We  first  look  at  the  special  case  when  u  —  w,  in  which  case  condition 
(10.9)  implies 


lim 

r—¥  00 


crMCr) 

rjrau(r]r) 


=  1. 


(10.10) 


We  refer  to  the  proof  of  Lemma  9.1,  where  we  introduced  the  asymptotically  equiv¬ 
alent  representation  (9.3)  of  the  slowly  varying  function  u.  Substituting  (9.3)  in 
(10.10),  we  see  that 


lim  exp 

r— koo 


a  —  e(t) 
t 


a  —  e(t) 
~~t 


=  1, 


or,  equivalently, 


a  —  e(t)  , 
lim  /  - -  dt 

r->°°  Jc,T  t 


=  0. 


(10.11) 


Pick  6  in  (0,  a).  Because  lim  Cr  =  lim  rjr  =  00,  there  exists  rg  such  that  for  r  >  rg 

r— >00  r— >00 

we  have  |e(f)|  <  8  whenever  t  >  min(<}r, r}T).  Thus, 


(a -5) 


Inf 

HO 

HO 

Ul 

1 

$ 

£ _ 

V 

Cr 

\J  Cr  t 

r  >  rs, 
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and  combining  this  last  inequality  with  (10.11)  we  obtain  the  desired  conclusion 

lim  rir/(r  =  1. 
r— >  oo 

In  general,  when  u  and  w  are  not  necessarily  equal,  we  note  the  easy  relation 


Cr  °u(Cr)  =  r(r  au{(r)  .  w{qr) 

T]rau(T]r)  rr]raw(r]r)  u(Tjr )  ’ 


(10.12) 


Condition  (10.9)  and  the  asymptotic  equivalence  of  u  and  w  together  imply  that 
the  relation  (10.10)  still  holds,  and  the  conclusion  Cr  ~  Vr  (r  —• ►  oo)  follows  from  the 
first  part  of  the  proof.  ■ 


Proposition  10.3  Consider  a  slowly  varying  function  u  :  1R+  — >  1R+,  and  let  a  >  1. 
For  any  sequence  {Cr,  r  =  1,2, . . .}  with  ^liinCr  =  oo  such  that  (10.9)  holds,  we  have 

( r  ~  r»w(r)  (r  — >  oo)  (10.13) 

for  some  slowly  varying  function  w  :  JR+  — >  1R._(- . 


Proof.  We  go  back  to  the  proof  of  Lemma  9.1,  where  we  introduced  the  asymp¬ 
totically  equivalent  representation  (9.3)  of  the  slowly  varying  function  u.  In  view  of 
Lemma  10.2,  it  suffices  to  consider  a  sequence  {Cr,  r  =  1,2, . . .)  determined  by  the 
relations 

rir-°-cexJj)'eMdt 


)  =  K.  r  >  r,  (10.14) 


with  r*  large  enough,  with  constants  A  >  0  and  c  >  0,  and  Borel  mapping  e  :  1R+  -> 
1R.  such  that  lini£(t)  =  0.  By  simple  manipulations,  we  can  write  (10.14)  in  the 

t— \oo 


equivalent  form 


Bra 


i  -  m 

t 


r  >  r* 


(10.15) 


with  i 

B^{jfAc)a  and  =  t~°'  (10.16) 

Hence,  by  Proposition  10.1,  for  large  enough  r  we  see  that  Cr  is  the  unique  solution 
x(y)  of  the  equation  (10.1)  with  y  =  Bra  .  By  the  second  part  of  Proposition  10.1, 
we  have 

Cr  =  x(Bra)  ~  Bray(Bra)  (r  — >■  oo),  (10.17) 
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and  the  desired  conclusion  is  now  immediate  once  we  note  that  the  mapping  w  : 
x  — »  Bv(Bxi  )  is  slowly  varying  whenever  v  is.  ■ 
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